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Relativistic Faddeev equations for three-body scattering at arbitrary energies are formulated in 
momentum space and in first order in the two-body transition-operator directly solved in terms 
of momentum vectors without employing a partial wave decomposition. Relativistic invariance is 
incorporated within the framework of Poincare invariant quantum mechanics, and presented in some 
detail. Based on a Malfliet-Tjon type interaction, observables for elastic and break-up scattering are 
calculated up to projectile energies of 1 GeV. The influence of kinematic and dynamic relativistic 
effects on those observables is systematically studied. Approximations to the two-body interaction 
embedded in the three-particle space are compared to the exact treatment. 

PACS numbers: 21.45+v 



I. INTRODUCTION 

Light nuclei can be accurately modeled as systems of nucleons interacting via effective two and three-body forces 
motivated e.g. by meson exchange. This picture is expected to break down at a higher energy scale where the physics 
is more efficiently described in terms of sub-nuclear degrees of freedom. One important question in nuclear physics 
is to understand the limitations of models of nuclei as systems of interacting nucleons. Few-body methods have been 
an essential tool for determining model Hamiltonians that describe low-energy nuclear physics. Few-body methods 
also provide a potentially useful framework for testing the limitations of models of nuclei as few nucleon systems, 
however this requires extending the few-body models and calculations to higher energy scales. There are a number 
of challenges that must be overcome to extend these calculations to higher energies. These include replacing the 
non-relativistic theory by a relativistic theory, limitations imposed by interactions fit to elastic scattering data, new 
degrees of freedom that appear above the pion production threshold, as well as numerical problems related to the 
proliferation of partial wave at high energies. In this paper we address some of these questions. We demonstrate that 
it is possible to now perform relativistic three-body scattering calculations at energies up to 1 GeV laboratory kinetic 
energy. The key elements of our success is the use of direct integration methods that avoid the use of partial waves 
and new techniques for treating functions of non-commuting operators that appear in the relativistic nucleon-nucleon 
interactions. 

During the last two decades calculations of nucleon-deuteron scattering based on momentum-space Faddeev equa- 
tions [l[ experienced large improvements and refinements. It is fair to state that below about 200 MeV projectile energy 
the momentum-space Faddeev equations for thrce-nuclcon (3N) scattering can now be solved with high accuracy for 
realistic two- and three-nucleon interactions. A summary of these achievements can be found in Refs. 0, @, ULIU- The 
approach described there is based on using angular momentum eigenstates for the two- and three-body systems. This 
partial wave decomposition replaces the continuous angle variables by discrete orbital angular momentum quantum 
numbers, and thus reduces the number of continuous variables to be discretized in a numerical treatment. For low pro- 
jectile energies the procedure of considering orbital angular momentum components appears physically justified due 
to arguments related to the centrifugal barrier and the short range of the nuclear force. If one considers three-nucleon 
scattering at a few hundred MeV projectile kinetic energy, the number of partial waves needed to achieve convergence 
proliferates, and limitations with respect to computational feasibility and accuracy are reached. It appears therefore 
natural to avoid a partial wave representation completely and work directly with vector variables. This is common 
practice in bound state calculations of few-nucleon systems based on variational Q and Green's function Monte Carlo 
(GFMC) methods @, 1, 1 G3 and was for the first time applied in momentum space Faddeev calculations for bound 
states in [ll| and for scattering at intermediate energies in Ref. fl2| . 

The key advantage of a formulation of the Faddeev equations in terms of vector variables lies in its applicability 
at higher energies, where special relativity is expected to become relevant. Poincare invariance is an exact symmetry 
that should be satisfied by all calculations, however in practice consistent relativistic calculations are more numerically 
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intensive, thus making their nonrelativistic counterpart a preferred choice. Furthermore, estimates of relativistic effects 
have been quantitatively small for 3N scattering below 200 MeV [TH, Q, EH with the exception of some breakup cross 
sections in certain phase space regions [l6| , indicating that at those energies non-relativistic calculations have sufficient 
precision. This is in part because in either a relativistic or non-relativistic model the interactions are designed to 
fit the same invariant differential cross section [l7| . which can be evaluated in any frame using standard kinematic 
Lorentz transformations, so model calculations are designed so that there are no "relativistic corrections" at the 
two-body level. Three-body interactions can even be chosen so the non-relativistic calculations fit both the two- and 
three-body invariant cross sections. This can be done in one frame and the invariance of the cross section fixes it in 
all other frames using standard relativistic kinematics. This, procedure has internal inconsistencies which show up if 
these models are used as input in larger systems, but they clearly indicate that the problem of identifying relativistic 
effects is more subtle than simply computing non-relativistic limits. In this paper we focus on differences between 
relativistic and non-relativistic calculations with two-body input that have the same cross section and use the same 
two-body wave functions [H, [H, HH ■ 

There are two primary approaches for modeling relativistic few-body problems. One treats Poincare invariance as a 
symmetry of a quantum theory, the other is based on quasipotential reductions [2l| of formal relations [12, H3| between 
covariant amplitudes. One specific realization of this approach is the covariant spectator approach of Ref. [2~4| . In 
this paper the relativistic three-body problem is formulated within the framework of Poincare invariant quantum 
mechanics. It has the advantage that the framework is valid for any number of particles and the dynamical equations 
have the same number of variables as the corresponding non-relativistic equations. Poincare invariance is an exact 
symmetry that is realized by a unitary representation of the Poincare group on a three-particle Hilbert space. The 
dynamics is generated by a Hamiltonian. This feature is shared with the Galilean invariant formulation of non- 
relativistic quantum mechanics. The Hamiltonian of the corresponding relativistic model differs in how the two-body 
interactions are embedded in the three-body center of momentum Hamiltonian (mass operator). The equations we 
use to describe the relativistic few-problem have the same operator form as the nonrelativistic ones, however the 
ingredients are different. 

In this article we want to concentrate on the leading order term of the Faddeev multiple scattering series within 
the framework of Poincare invariant quantum mechanics. The first order term contains already most relativistic 
ingredients which, together with the relativistic free three-body resolvent, gives the kernel of the integral equation. 
We want to understand essential differences between a relativistic and nonrelativistic approach already on the basis of 
the first order term. As simplification we consider three-body scattering with spin-independent interactions. This is 
mathematically equivalent to three-boson scattering. The interactions employed are of Yukawa type, and no separable 
expansions arc employed. In order to obtain a valid estimate of the size of relativistic effects, it is important that the 
interactions employed in the nonrelativistic and relativistic calculations are phase-shift equivalent. To achieve this 
we employ here the approach suggested by Kamada-Glockle [l8j . which uses a unitary rescaling of the momentum 
variables to change the nonrelativistic kinetic energy into the relativistic kinetic energy. 

This article is organized as follows. Section II discusses the formulation of Poincare invariant quantum mechanics, 
and section III discusses the structure of the dynamical two- and three-body mass operators. Scattering theory formu- 
lated in terms of mass operators is discussed in section IV. The formulation of the Faddeev equations and techniques 
for computing the Faddeev kernel are discussed in section V. Details on kinematical aspects of how to construct the 
cross sections is given in section VI. In Sections VIII and IX we present calculations for elastic and breakup processes 
in the intermediate energy regime from 0.2 to 1 GeV. Our focus here is to compare different approximations to the 
embedded interaction with respect to the exact calculation. Our conclusions are in Section X. Two Appendices are 
devoted to relating the transition matrix elements based on mass operators to the invariant amplitudes with the 
conventions used in the particle data book and expressing the invariant cross section and differential cross sections 
worked out directly in laboratory-frame variables. 

II. POINCARE INVARIANT QUANTUM MECHANICS 

Symmetry under change of inertial coordinate system is the fundamental symmetry of Poincare invariant quan- 
tum mechanics. In special relativity different inertial coordinate systems are related by the subgroup of Poincare 
transformations continuously connected to the identity. In this paper the Poincare group refers to this subgroup, 
which excludes the discrete transformations of space reflection and time reversal. Wigncr [25| proved that a necessary 
and sufficient condition for quantum probabilities to be invariant under change of inertial coordinate system is the 
existence of a unitary representation, U(A, a), of the Poincare group on the model Hilbert space. Equivalent vectors, 
l^) and in different inertial coordinate systems are related by: 



(2.1) 
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In Poincare invariant quantum mechanics the dynamics is generated by the time evolution subgroup of U(A,a). 
The fundamental dynamical problem is to decompose U(A, a) into a direct integral of irreducible representations. 
This is the analog of diagonalizing the Hamiltonian or time-evolution operator in non-relativistic quantum mechanics. 
The problem of formulating the dynamics is to construct the dynamical representation U(A, a) of the Poincare group 
by introducing realistic interactions in the tensor products of single particle irreducible representations in a manner 
that preserves the group representation property and essential aspects of cluster separability. The solution to this 
non-linear problem is achieved by adding suitable interactions to the Casimir operators of non-interacting irreducible 
representations of the Poincare group. 

Since irreducible representations of the Poincare group play a central role in both the formulation and solution of 
the dynamical model, we give a brief summary of the construction of the irreducible representations that we use in 
this paper. The Poincare group is a ten parameter group that is the semidirect product of the Lorentz group and 
the group of spacetime translations. Spacetime translations are generated by the four momentum operator, P' 1 , and 
Lorentz transformations are generated by the antisymmetric angular momentum tensor, J^ v . 

The Pauli Lubanski vector is the four vector operator defined by 

W» = ~e» a ^P a Jp y . (2.2) 

The Casimir operators for the Poincare group are 

M 2 = -rj^P^P" = H 2 P P (2.3) 

and 

W 2 = -rj^W^W = M 2 j 2 (2.4) 

where rj^ is the Minkowski metric, M is the mass operator, H is the Hamiltonian, P is the linear momentum, and 
j 2 is the spin. 

Positive-mass positive-energy irreducible representations are labeled by eigenvalues of the mass M and spin j . 
Vectors in an irreducible subspace are square integrablc functions of the eigenvalues of a complete set of commuting 
Hcrmitian operator- valued functions of the generators -P M and J MI/ . In addition to the two invariant Casimir operators, 
it is possible to find four additional commuting Hcrmitian functions of the generators. For each of these four commuting 
observables it is possible to find conjugate operators. These conjugate operators, along with the eigenvalues of 
the Casimir operators, fix the eigenvalue spectrum of the four commuting Hcrmitian operators. The irreducible 
representation space is the space of square integrable functions of the eigenvalues of the four commuting operators. 
The generators can be expressed as functions of these four operators, their conjugates, and the Casimir invariants 

[IS EE HI. 

In this paper we choose the four commuting operators to be the three components of the linear momentum and 
the z component of the canonical spin operator. In this representation the four conjugate operators are taken as the 
partial derivatives of the momentum holding the canonical spin constant (Newton- Wigner position [2^] operator) and 
the x component of the canonical spin. While j x is not exactly conjugate to j z , the two operators generate the full 
SU(2) spin algebra. The corresponding eigenstates have the form 

|P,M> : = l(™»i)P>A*)- ( 2 -5) 

The mass m spin j irreducible representation of the Poincare group in this basis is determined from the group 
representation property and the action of rotations, spacetime translations and canonical boosts on the zero momentum 
eigenstates: 

U(R,0)\O^) = \O^')Dl^(R) (2.6) 
U(l,a)\0,fi} = e- ia ° m \0,fi) (2.7) 



( 2 - 
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where in these equations R is a rotation, D 3 , (R) is the standard 2j + 1 dimensional unitary representation of SU(2), 
a = (a , a) is a displacement four vector, B{p m ) is the rotationless Lorentz boost (canonical boost) that transforms 
(to, 0) to p m ■= (E Pm ,p), 

(%))V=(^" I+ »t ), (2.9) 

\ 1 ' m(m+E Pm ) / 



and E Pm = y/ m 2 + p 2 . That the magnetic quantum number remains invariant in (|2.8I) under the rotationless boost 
(|2.9p is the defining property of the canonical spin. 

The multiplicative factor on the right side of (|2.8p is fixed up to a phase by unitarity and the normalization 
convention 

(p',fM'\p, f x) = S(p'-p)d^ fl . (2.10) 
With these choices the action of an arbitrary Poincare transformation on these states is given by 




W mj (A,a)|p, M ) = \p'^') x -^D 3 (R w {A,p m ))e^ a , (2.11) 



where R w (A,p m ) is the (standard) Wigner rotation, 

R w (A, Pm ) :=B- 1 (Ap m )AS(p m ), (2.12) 

and p' m = Ap m . Since each of the elementary transformations (|2.6I2.7|2.8|) are unitary, it follows that (|2.11|) is unitary. 
Since every basis vector can be generated from the [i = j and P = basis vector using equations (|2.6I2.7I2.8[) . 
representation (|2.1ip is also irreducible. 

The mass m spin j irreducible representations that are used in this paper have the form (|2.11| . The mass m spin 
j irreducible representation space with basis (|2.5p is denoted by H m j . 

The Hilbert space for the-three nucleoli problem is the tensor product of three one-nucelon irreducible representation 
spaces: 

n = n mh ®n mh ®n mh . (2.13) 

In this paper all nucleons are assumed to have the same mass, to. 

The non-interacting unitary representation of the Poincare group on Ti is the tensor product of three one-nucleon 
irreducible representations: 

Wo(A,o) =U mh (A,a)®U m i(A,a)®U m i(A,a). (2.14) 

As in the case of rotations, the tensor product of irreducible representations of the Poincare group is reducible. 
The tensor product of three irreducible representations can be decomposed into a direct integral of irreducible repre- 
sentation using Clebsch-Gordan coefficients for the Poincare group. The Clebsch-Gordan coefficients for the Poincare 
group are known [28|, HH, [32j . As in the case of rotations, the Poincare Clebsch-Gordan coefficients are basis dependent 
and the three-body irreducible representations can be generated by pairwise coupling. The Poincare Clebsch-Gordan 
coefficients can be computed by using SU(2) Clebsch-Gordan coefficients to decompose three-body zero momen- 
tum eigenstates states into irreducible SU(2) representations. Three-body irreducible basis vectors are generated by 
applying Eqs. ()2.6I2.7I2.8[) to the zero-momentum SU{2) irreducible representations. 

The resulting irreducible three-body basis depends on the order of the coupling. In the basis of eigenstates of 
the three-body linear momentum and canonical spin the irreducible eigenstates are labeled by eigenvalues W of the 
three-body invariant mass, Mq, the three-body canonical spin, j (for simplicity of notation we use the same label 
j for the total canonical spin of the three-body system and the single particle canonical spin), the total three-body 
momentum, P, the z component of the three-body canonical spin, /x, and invariant degeneracy quantum numbers, d, 
that distinguish multiple copies of the same irreducible representation: 

|(W,j),P, M ,d). (2.15) 

For two-particle systems the degeneracy quantum numbers d are discrete (for example they may be taken to be 
invariant spin and orbital angular momentum quantum numbers) while for more than two particles the degeneracy 
quantum numbers will normally include invariant sub-energies, which have a continuous eigenvalue spectrum. In 
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addition to the appearance of the degeneracy quantum numbers, the eigenvalue spectrum of the free invariant mass 
operator, M , is continuous. 

These states transform as mass W spin j irreducible representations of the Poincare group under Uq(A, a): 



U (A, o)| (W, j), P, n, d) = |(W, j), P fj, d)J ^D 3 ^ (i?»(A, fh)) e tP «' a , (2.16) 



where 



= AP W P w := (v/W2 + P2,P) = (Ep„,P). (2.17) 



The quantities W, j, d are invariants of the representation (|2.16|) of Uo(A, a). 

Because the Poincare group allows time evolution to be expressed in terms of spatial translations and Lorentz 
boosts, when particles interact, consistency of the initial value problem requires that the unitary representation of 
the Poincare group depends non-trivially on the interactions. The construction of U(A, a) for dynamical models is 
motivated by the example of Galilean invariant quantum mechanics. The non-relativistic three-body Hamiltonian has 
the form 

H =& g +k (2 ' 18) 

where the Casimir Hamiltonian, h, is the Galilean invariant part of the Hamiltonian and M g is the Galilean mass. In 
the non-relativistic case interactions are added to the non-interacting Casimir Hamiltonian ho : 

h = h + V nr (2.19) 

where the Galilean invariance of h requires that the interaction V nr is rotationally invariant and commutes with and 
is independent of the linear momentum P. This means that in the corresponding non-relativistic basis 

(h,p 1 j )M ,d|y„ r |h^p^j^/i^ f ^o=^(p-po%'V( h ^ll^rllh',rf')■ (2-20) 

where h is an eigenvalue of ho- 

In the Poincare invariant case Eq. (|2. 18[) is replaced by 



H = VP 2 + M 2 (2.21) 

where M in Eq. (|2.21[) plays the same role as the Casimir Hamiltonian h in Eq. (|2.18|) . The corresponding free 
relativistic Hamiltonian is Ho = \fV 2 + Mq. In what follows A denotes the eigenvalue of M to distinguish it from 
the eigenvalue W of Mo. 

A Poincare invariant dynamics can be constructed by adding an interaction to the non-interacting Mo that commutes 
with and is independent of P and j z : 

M = M + V r . (2.22) 
In the non-interacting irreducible basis (|2.15[) these conditions require interactions of the form 

((w,^,p,M,rf|K|(w^/),p^ / i^do = ( 5(p-po%'V(w,d!lK• ? !lw',rf')• (2.23) 

The dynamical problem is to find simultaneous eigenstates of the commuting operators M, P, j z , and j 2 . This is 
done by diagonalizing M in the irreducible basis (|2.15j) . The cigenfunctions of M have the form 

((W, j), P, fx, d\(X, /), P ', //) = *(P - P ')%< V (W, d, j\X,j) (2.24) 

where the eigenfunctions (W, d,j\X, j) are solutions of 

\N(\N,j,d\X,j} + J2 I (W,d||Vv'||W',rfVW , ^'<W',. ? ,d'|A,j) - \(\N,j,d\X,j), (2.25) 

with mass eigenvalue A. For two-particle systems the degeneracy quantum numbers d arc discrete (for example they 
may taken to be invariant spin and orbital angular momentum quantum numbers) while for more than two particles 
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the degeneracy quantum numbers will normally include invariant sub-energies, which have a continuous eigenvalue 
spectrum. 

Because {Mo, j 2 , j z , P, j x , — iVp} have the same commutations relations as I M, j 2 , j z , P, j x , — iVp}, if the dynami- 
cal Poincare generators are defined as the same functions of these operators [27|,[28[ with Mq replaced by M, it follows 
that the simultaneous eigenstates 

l(A,j),P,M) (2-26) 

of M, j 2 ,P,j z transform as a mass A spin j irreducible representation of the Poincare group. 

Since these eigenstates are complete, this defines the dynamical representation of the Poincare group on a basis by 

W(A,a)|(A,j),P,M) = |(A,j),P', A 0^/|^^(^(A,PA))e^- a , (2.27) 

where 

E Px = \/A 2 + P 2 , Pa := {E P ,P), P' x = AP A . (2.28) 

This shows that the solution to the eigenvalue problem (|2.25[) provides the desired decomposition of the dynamical 
unitary representation of the Poincare group into a direct integral of irreducible representations. 

The appearance of the mass eigenvalue A on the right side of Eq. (|2.27[) indicates the interaction dependence of this 
representation. It can happen, for a given choice of irreducible basis, that the coefficient on the right hand side of (|2.27p 
is independent of the mass eigenvalue for a subgroup of the Poincare group. For the basis (|2.15[) . of linear momentum 
and canonical spin eigenstates, both translations and rotations have this property. These transformations generate a 
three-dimensional Euclidean subgroup of the Poincare group that is independent of the interaction. An interaction- 
independent subgroup is called a kinematic subgroup; the three-dimensional Euclidean group is the kinematic subgroup 
is associated with an instant-form dynamics [301 ] . 



III. MASS OPERATORS 



Next we discuss the structure of mass operators for the two and three-body problems. We pay particular attention 
to issues related to representations of these operators that are suitable for computations without using partial waves. 

The construction of the dynamics in (|2 . 25[) and (|2.27[) adds an interaction to the mass Casimir operator of a non- 
interacting irreducible representation of the Poincare group to construct an interacting irreducible representation. The 
role of the spin in the structure of the irreducible representations suggests that this construction requires a partial 
wave decomposition, however partial waves are not used in the calculations that follow. 

The spin in the relativistic case is obtained by coupling the single particle spins and orbital angular momentum 
vectors. While the form of the coupling is more complex than it is in the non-relativistic case, the final step involves 
coupling redefined spins and orbital angular momenta with ordinary SU(2) Clebsch- Gordon coefficients. Undoing 
this coupling leads to a representation of the dynamical operators that can be used in a calculation based on vector 
variables. 

The first step is to construct redefined vector variables that can be coupled to obtain the spin. To understand the 
transformation properties of these operators note that the magnetic quantum number in equation (|2.27p undergoes a 
Wigncr rotation when the system is Lorentz transformed. If the spin of the representation is obtained by coupling the 
redefined particle spins and orbital angular momenta with SU(2) Clcbsch-Gordan coefficients, then all of the spins 
and relative angular momenta must also transform with the same Wigner rotation. 

To illustrate how to construct momentum operators that Wigner rotate under kinematic Lorentz boosts consider a 
pair of noninteracting spinless particles. The total four momentum Pm of this system is the sum of the single particle 
four momenta 

Pm = p mi +Pm 2 - (3.1) 

Define the operator k by 

k:=B- 1 {P Mo )p m ^ (3.2) 

where B~ 1 (P]\i a ) is interpreted as a 4 x 4 matrix of operators. If both Pm and p mi are transformed with a Lorentz 
transformation A, then k rotates with the same Wigner rotation that appears in Eq. (|2.16|) 



k' = B-^P^M = B-\AP Mo ) A Pl = j B- 1 (AP Mo )AB(P Mo )B- 1 (Pmo)pi = ^(A, P Mo ) k. (3.3) 
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It is due to the operator nature of B~ 1 {Pm ) that k docs not transform as a 4-vector. 

The tensor product of single particle basis vectors |pi, P2) can be replaced by a basis |P, k) using a variable change. 
In this basis Eq. (|3.3jl implies 



Wo(A,o)|P J k) = |P',i?4A,/V)k)J^e iP w- a P^ = AfW. (3.4) 

V h Pw 



where k is the eigenvector of the space components of the operator (|3.2[) . This shows that k undergoes the same 
Wigncr-rotation as the two-body canonical spin. 

If the two particles have spin, this single particle spins need to be Wigner rotated before they can be coupled [28j . 
The spins obtained this way are called constituent spins. The constituent spins undergo the same Winger rotations 
as k but they are not 1-body operators and do not have natural couplings to the electromagnetic interaction. In this 
paper we only consider spinless interactions. In this case the constituent spins can be ignored. Their only effect is 
that they impact the permutation symmetry of orbital- isospin part of the wave functions. 

The magnitude of k is an invariant that fixes the two-body invariant mass eigenvalue W: 

W = 2y/k 2 +m 2 . (3.5) 
If the vector part of k is expanded in partial waves 

|P,|kU M ) := J |P,k)y^(k)dk (3.6) 

then 

U (A, a)|P, |k|, j, n) := £ |P ', |k|, j, //) ^^^(^(A, Pw))e iP ™- a (3.7) 

transforms like (|2. 1 1|) . 

In the representation |P,k) Eq. (|2.23|) is satisfied for a spinless interaction of the form 

(P ! k|K|P',k')=5(P-P')(k||y r ||k') (3.8) 

with a rotationally invariant kernel 

(i?k||Vv||i?k') = (k||Vv||k'). (3.9) 

If the interaction includes nucleon spins, the rotational invariance must be generalized to include rotationally invariant 
contributions involving the constituent spins. 

Next we consider the three-body problem, where Wo (A, a), Mo and W are now associated with the three nucleon 
system. In the three-body system vector operators that Wigner rotate are the Poincare- Jacobi momenta and three- 
body constituent spins. The Jacobi momenta are obtained from the non-relativistic Jacobi momenta by replacing 
Galilean boosts to the zero momentum frame of a system or subsystem by the corresponding non-interacting Lorentz 
boosts. In these expression the boosts are considered to be matrices of operators. The replacements are 

ky EE B; 1 (P)( Pl - p,) = Pi - Pj - (P '+ Pj) K - mj ) — ► EE B-\P mQ )&* - Pj ) (3.10) 

y (rrii + mj) 

q t EE S- 1 (P) Pi = Pi - ^-rm — ► q, ee B-\P Mo ) Pi . (3.11) 

Mg 

where mo,y is the invariant mass of the noninteracting two-particle system. 

The only complication in the three-body case is that when the single particle momenta undergo Lorentz transfor- 
mations the variables ky , q 4 experience different Wigner rotations 

q, -* R w (A,P Mo H (3-12) 
kjj ► R w (A, P ma ij )kjj . (3.13) 

Because of the different Wigner rotations, the angular momenta associated with q^ and ky cannot be consistently 
coupled with SU (2) Clebsch- Gordon coefficients. To fix this redefine ky — > ky by replacing all of the p^s in (|3.10D 
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by the corresponding t^s. Then when the single particle variables are Lorentz transformed, the all Wigner rotate 
with a rotation R' . This means that the redefined k^ transform as R W {R' , qij), where q^ = + qj. But the defining 
property of the canonical boost (|2.9|) is R W (R', qij) — R' which means that both and ky undergo the same Wigner 
rotation, i?', when the single-particle variables are Lorentz transformed. 

Only two of the six vector variables, q; and ky, are linearly independent. Any two of these variables along with 
P can be used to label three-body basis vectors. Following non-relativistic usage, the single particle momenta are 
replaced by the independent variables 

{P,q*,k«} (3-14) 
where k ^ The single particle basis vectors are replaced by 

lr> l \ I \ i <5(Pli P2j P3) 1I/2 ^oicn 

|P,q fe ,ky) = | Pl ,p 2 ,p 3 )|_ —| / . (3.15) 

1 Hk, K-ij ) 

The above definitions imply the desired transformation property 

Wb(A,a)|P,q*,ky) = |P',i? lu (A,PW)q fc ,i?. lu (A,PW)k y ) J^e jp w- a (3.16) 

V 

where P' = AP. 

The operators q^ and k; are functions of the single particle momenta, and are thus defined on states with any total 
momenta, not just on three-body rest states. This is similar to the mass operator, which is also defined on states of 
any total momentum. 

Next we discuss the structure of the mass operators that will be used in the two and three-body problems. The 
mass operators mo,jj and Mq for the non-interacting two- and three-body systems can be expressed in terms of the 
operators ky, ky, and q^ as 



mo,ij = \jm\ + k£ + yjm) + k?. (3.17) 

and 



where mo,ij in (|3. 1 8[) replaces k^ by k^. 

When two-body interactions are added to rriQ.ij the interacting two-body mass operator becomes: 

my = m ,y + Vij (3.19) 

where in this basis (|2.23[) becomes 

(py ,kij\vij\t>ij ',kij ') = Sfaij - '){&4j\\vij\\kij ')• (3.20) 

Cluster properties determine how the two-body interactions enter the three-body mass operator. In order to obtain 
a three-body scattering operator that clusters properly, it is enough to replace kjj by ky in the two-body interaction 
and include the modified two-body interaction in the three-body mass operator as follows: 

M tj = y/imort+Virf + qfc + ^ m 2 k + <$ (3.21) 

where 

(P, q fc , kijlvij |P ', q fc ', kij ') = £(P - P ')*(q* - <1* ')( k «ll^'l|ky ') (3.22) 

and the functional form of the reduced kernel (ky||?jy||ky ') is identical in (|3.20p and (|3 . 22[) . and it must be a 
rotationally invariant function of its arguments k^ and kjj ' (resp. ky and kjj ') . 
The interacting three-body mass operator is then defined by 

M := M + Via + ^23 + V 31 (3.23) 
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where the two-body interactions embedded in the three-particle Hilbcrt space 32j are 



Vij := My - M = ^(mo^+v^y+ql - ^m^.+qf.. (3.24) 

The sum of the interactions is consistent with the constraint (|2.23|) since each of the two-body interactions Vij in (|3.23|) 
is consistent with (|2.23|) . The dynamical problem is to find eigenstates of M in the basis (|3.15|) . The technical challenge 
for the numerical solutions of the three body-problem is the computation of the embedded two-body interactions, Vij, 
which requires computing functions of the non-commuting operators mo.y and Vij. 

While this dynamical model leads to an S'-matrix that satisfies cluster properties, the constructed unitary repre- 
sentation of the Poincare group only clusters properly when P = 0. Since the ^-matrix is Poincare invariant, this 
is sufficient for computing all bound state and three-body scattering observables, however additional corrections are 
required if the three-body eigenstates are used to compute electromagnetic observables. 



IV. RELATIVISTIC SCATTERING THEORY 



Our interest in this paper is the computation of scattering cross sections for two-body elastic scattering and breakup 
reactions in Poincare invariant quantum mechanics. The formulation of the scattering theory using dynamical mass 
operators for the Poincare group is outlined below. For a more complete discussion see (28j . 

The multichannel scattering matrix is calculated using the standard formula 

Sa0 = (^(0)|^(0)) (4.1) 

where (3 is the initial deuteron-nucleon channel and a is either a deuteron-nucleon or three nucleon channel. 

The scattering states that appear in Eq. (|4.ip are defined to agree with states of non-interacting particles long 
before or long after the collision 

lim || e -^|*±(0))-e-^*|$±(0))||=0. (4.2) 

t — >±oc 

In this paper the ± on the scattering states and wave operators indicates the direction of the time limit 
(— =past/+ =future), which is opposite to the sign of it. 

In the breakup channel I^^O)) is a normalizable Hubert space vector. In the nucleon-deuteron channels (0)) 
has the form 

(P,q 4 ,k, fe |$±(0)) = <Mk, fc )/(q 4 ,P) (4.3) 

where (poO^-jk) is the deuteron wave function and /(q^P) is a unit normalized wavepacket describing the state of a 
free deuteron and nucleon at time zero. 

The asymptotic and interacting scattering states arc related by the multichannel wave operators 

|*±(0))=n a± (H > H a )|*±(0)) (4.4) 

where the multichannel wave operators are defined by the strong limits 

il a± = lim e zHt e - lHat (4.5) 

t— »±oo 

on channel states. The multichannel scattering operator can then be expressed in terms of the wave operators as 

S afj = nl + (H,H a )n _(H,H ). (4.6) 

In the three-body breakup channel, a = 0, 



H a = H Q = ^M$ + P2. (4.7) 
In channels, a = (ij), with an incoming or outgoing deuteron, 

H a =H ij =H + V Ha , (4.8) 

where 



V Ha = JMfj + P 2 - JM§ + P 2 (4.9) 
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is the interaction between the nucleons in the deuteron in the three-body Hamiltonian. We use the notation M a 
to denote M for the breakup channel or My for the nucleon-deuteron channel. In the non-relativistic case the 
Hamiltonian. which generates time evolution in the asymptotic conditions, is normally replaced by the Casimir 
Hamiltonian, h. This can be done because P 2 /2Af 9 appears linearly in both the interacting and non- interacting 
Hamiltonians and commutes with the interactions. In the relativistic case the interaction in the Hamiltonian is 
different from the interaction in the mass operator, and the kinetic energy enters the mass non-linearly. In the 
relativistic case the wave operators can still be expressed directly in terms of the mass operators. The justification 
is the Kato-Birman invariance principle [HI, |34| which implies that H and H a can be replaced by a large class of 
admissible functions of H and H a in the wave operators; specifically 

M = \/H 2 - P 2 (4.10) 

is in the class of admissible functions. This gives 

n a± = km e iMt e- iM ~ t = lim e^e" 4 ^* (4.11) 

t — >±oo t — >±oo 

which leads to an expression for the multichannel scattering operator [26| expressed directly in terms of the mass 
operators: 

S a p = lim e iM a r e -iM{r+r-) e iM s r^ (412) 
t,t' — >oc 



If these limits are computed in channel mass eigenstates |a) and \0) of M a and Mp the result is 

(a\S\fl = (a\0) - 2mS(\N a - W a )(a|T Q/3 (W Q + i0 + )\p), (4.13) 

where 

T al3 (z) = V^ + V a (z-M)- 1 V , (4.14) 

and 



V a = M - M a = M — Mij (4.15) 

for two cluster n — d channels and 

V a = M — M (4.16) 

for the breakup channel. Here W a and are the eigenvalues of M a and Mp in the channel eigenstates |a) and \(3). 
The first term in Eq. (|4.13[) is identically zero if the states \a) and \/3) correspond to different scattering channels. 

Compared to the standard expression that is based on using the Hamiltonian, in Eq. (|4.14|) the interactions are 
expressed as differences of mass operators rather than Hamiltonians, the resolvent of the Hamiltonian is replaced 
by the resolvent of the mass operator and the energy conserving delta function is replaced with an invariant mass 
conserving delta function. 

The translational invariance of the interaction (|4.15j) requires that 

(P, • • • \T^{z)\ ■ ■ ■ ,P'> = 50? - P')(- • • \\T afj (z)\\ ■ ■ •). (4.17) 

Given the momentum conserving delta function, the product of the momentum and mass conserving delta functions 
can be replaced by a four-momentum conserving delta function and a Jacobian: 

5(W Q - W^)<5(P Q - P p ) = S\P a - Pp) 



where E a = ^/W 2 + P 2 . 

The representation (|4.13p of the scattering matrix can be used to calculate the cross section. The relation between 
the scattering matrix and the cross section is standard and can be derived using standard methods, such as the ones 
used by Brenig and Haag in [35| . The only modification is that in the usual expression relating the cross section to 
the transition operator, the transition operator is the coefficient of — {2ir)i8(E a — Ep). Thus to compute the cross 



dE 



dM 



5 4 (P a - Pp) 



Wo 



(4.18) 
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section it is enough to use the standard relation between T and S with the channel transition operator being replaced 
by 

x(...||r^(z)||...) (4.19) 
The resulting expression for the differential cross section for elastic scattering is given by : 

da = ^\(p d ,p n \\T a f :i \\p , d ,p' n )\ 2 ^5 i ( Pd +p n -p' d -p' n ) dp n dp d . (4.20) 

V nd L Q 

and for breakup reactions the formula is replaced by 

da = I <Pi> P2- Pall^llPd, P'n) T % S 4 ( Pl +P2+ P3 - P ' d - p' n ) d Pl dp 2 dp 3 . (4.21) 

V nd 6 - L ct 

These relations are normally given in terms of single particle momenta while the transition matrix elements are 
evaluated in terms of the Poincare-Jacobi momenta. The transformation relating these representations involves some 
Jacobians. These are discussed in the sections on calculations. 

Except for the factor W^/E^, Eq. (|4.20|) is identical to the corresponding non-relativistic expression. The additional 
factor of W^/E^ arises because we have chosen to calculate the transition operator using the mass operator instead 
of the Hamiltonian. The differences in these formulas with standard formulas are (1) that the transition operator is 
constructed from the difference of the mass operators with and without interactions and (2) the appearance of the 
additional factor of W^/E^ which corrects for the modified transition operator. This factor becomes 1 when P = 0. 

The differential cross section da is invariant. Equations (|4.20[) and (|4.2ip can be expressed in a manifestly invariant 
form. The relation to the standard expression of the invariant cross section using conventions of the particle data 
book [36] is derived in the appendix where we also outline the proof of Eq. (|4. 1 3|) . 



V. INTEGRAL EQUATIONS 

The dynamical problem is to compute the three-body transition operators T a ^{z) that appear in Eqs. (|4.20p and 
(I4.21[) and use these to calculate the cross sections. It is useful to replace the transition operators (|4. 14[) by the 
on-shell equivalent AGS [38!] transition operators: 

U af3 (z) := 5 afi {z - M a ) + T a0 (z). (5.1) 

When z is put on the energy shell and evaluated on the channel eigenstate for the channel a the first term of (|5.1[) 
vanishes. The AGS operators are solutions of the integral equation 

U a/3 (z) := S a/3 (z - M ) + J2 ~ M 7 )- 1 C/'" 9 (z), (5.2) 

7 

where V 1 := M 1 — Mq = Vij are the embedded two-body interactions given in (|3.24p . and the sum is over the three 
two cluster configurations, 1 = (1; 23), 2 = (2; 31), 3 = (3; 12), with each cluster labeled by the index 7. 

When the particles are identical this coupled system can be replaced by an equation for a single amplitude, 

U(z) = P(z - M ) + PVl(z - A/ 1 )^ 1 C/(z), (5.3) 

where we chose without loss of generality to single out the configuration (1;23). In this case the permutation operator 
P is given by P = P12-P23 + Pi3p23- This solution can be used to generate the breakup amplitude 

U = (I + P)V 1 (z-M 1 )- 1 U(z). (5.4) 

The AGS operators U(z) and Uq(z) can be expressed in terms of the solution T(z) of the symmetrized Faddeev 
equations 

U(z) = P(z - M ) + PT(z) 

Uo(z) = (l + P)T(z), (5.5) 
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where T(z) is the solution to 

T(z) = Ti(z)P + T l {z)P{z - M )- 1 T(z). (5.6) 

where the operator T\(z) is the two-body transition operator embedded in the three-particle Hilbert space defined as 
the solution to 

T 1 {z) = V 1 +V 1 (z-M )- 1 T 1 {z), (5.7) 

where V\ = V23 = M23 — Afo- 

The first order calculation that we perform in this paper is defined by keeping the first term of the multiple scattering 
series generated by (|5.6|) : 

T{z)^T 1 {z)P (5.8) 



U(z)^P(z-M Q ) + PT 1 (z)P U (z)^(l + P)T 1 {z)P, (5.9) 

Because the embedded interactions V 1 ()3.24|1 in the AGS equations arc operator valued functions of the non- 
commuting operators TOo,ij and uy, their computation, given Vij as input, presents additional computational challenges 
To compute the kernel note that it can be expressed as 

Vy(z - M 7 ) _1 = T 7 (z){z - Ma)' 1 (5.10) 

where 

T 1 {z) = V 1 + V 1 {z-M 1 )- 1 V 1 (5.11) 

In this paper we compute this kernel using a method that exploits the relation between the two-body transition 
operator and the operator (|5. 1 1|) . Because My and rriij have the same eigenvectors it follows that [2(3] 

www - M.)-ip,„ k ) - «tp' - w,- - <» (Si^tra) ( " 2) 

where 

z' = M (q,k')+iO + 5' = m 0ij (k') + i0+. (5.13) 

In the AGS equation, this kernel is needed for all values of z, while equation (|5. 12[) only holds for z = z' . The kernel 
for an arbitrary z can be computed by using the first resolvent equation which leads to integral equation 

T 7 (z) = T y (z') + T 1 (z)G (z)(z' - z)G (z')T y (z') (5.14) 

where G(z) = (z — Mq)^ 1 which can be used to calculate T 7 (z) from T 7 (z') for all z 7^ z'. 



VI. RELATIVISTIC FORMULATION OF THREE-BODY SCATTERING 



In the scattering of three particles interacting with spin independent interactions, there are two global observables, 
the total cross section for elastic scattering, a e i, and the total cross section for breakup, o\, T . These can be computed 
using fOU]) and flgTlTJ) . In this section we discuss the kinematic relations needed to compute these quantities in more 
detail. 

If we replace the transition operators by the corresponding symmetrized AGS operators, use the identities 

n d = J |P,q, </3 d )dPdq(P,q, tp d \ = /"|p<i,Pn,Vd)^Pd^Pn(Pd,Pn,yd| (6.1) 



and 



I = J |P,k,q)dP(ikdq(P,k,q| = J |pi, p 2 , p 3 )dpidp 2 dp 3 (pi, P2, P3I, (6.2) 
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and evaluate the initial state and v b t in the center of momentum frame, (|4.20[) becomes 

(2^) 4 



Pel 



I'bt 



for elastic scattering and 



3 v bt 



dq.6QNf-\Ni)\(tp d ,q\\U\\(p d ,qo) 



dq dk<5(W / -Wi)|(k,q||C/ ||^ 1 qo) 



(6.3) 



(6.4) 



for breakup. 

Here Wi(W/) are the invariant masses eigenvalues of the initial (final) state, qo is the Poincare-Jacobi momentum 
between the projectile and the target, and k and q the Poincare-Jacobi momenta for a given pair and spectator defined 
in the previous section. The permutation operator in (|5.3[) only includes three of the six permutations of the three 
particles; the other three independent permutations are related by an additional transposition that interchanges the 
constituents of the deuteron, which is already symmetrized. This accounts for replacement of the statistical factor 
1/3! in (|4~20l) by the factor of 1/3 in (fO]) . 

Using relativistic kinematical relations the integral over |q| in Eq (|6.3p can be done using the invariant mass 
conserving delta functions with the result 



(7,1 



(2*0 



(6.5) 



The quantities W and qo are determined by the laboratory kinetic energy Ei ab of the incident nucleon. First note 

W 2 = (to + m d f + 2m d E lab . (6.6) 

The nucleon rest mass is given by to, the rest mass of the deuteron is m d = 2m — e dl where e d is the deuteron binding 
energy. The Poincarc Jacobi momentum between projectile and target, qo, is related to E\ ab by 



m d Ei ab 
W 2 



(E lab + 2m) 



(6.7) 



In the nonrelativistic case the phase space factor under the integral of Eq. (|6.5j) reduces to (2to/3) 2 . 

It is also necessary to compute the transition matrix elements that appear in (|6.3[) and (|6.4[) . The momenta of 
the three particles can be labeled either by single-particle momenta pi, P2, and p 3 , or the total momentum P and 
the relativistic Poincare-Jacobi momenta of Eqs. (|3. 10|) and (|3.11|) with the p ?; replaced by qi. The explicit relations 
between the three-body Poincare Jacobi momenta are 



-(qj + Qk) 



E, - E k 



k = ki = k jk = 1 (qj - q k ) - - (qj + q k ) ( 

2 2 + Ek + y/(Ej + Ek) 2 - (qj + q k ) 2 



(6.8) 



where Ei = E(qi) = y/ to 2 + qi 2 . For nonrelativistic kinematics the second term in k, being proportional to the total 
momentum of the pair (j,k), vanishes. In addition, the transformation from the single particle momenta pi to the 
Poincare-Jacobi momenta has a Jacobian given by 



|Pl,P2,P3) = 



9(P,k,q) 



<9(P1,P2,P3) 



1/2 



|P,k,q) 



where for P = the Jacobian becomes 



9(P,k,q) 



<9(P1,P2,P3) 



1/2 



|p=o 



V(E(q 2 ) + j^(q 3 )) 2 - q 2 (E(q 2 ) + E(q 3 )) ' 
AE(q 2 )E(q 3 ) 



1/2 



= 71 (q; q 2 q 3 ) 



(6.9) 



(6.10) 



In the above expression we chose without loss of generality particle 1 as spectator. The difference between the 
relativistic and non-relativistic Jacobi momenta in Eqs. (|3.10p and (|3.1ip are relevant for the calculation of the 
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permutation operator P in Eqs. (|5.5[) and (|5.6|) . The matrix elements of the permutation operator are then explicitly 
calculated as 



<k',q'|P|k,q) = N(q',q) 



«5(k' - q - lq' C(q, q'))<5(k + q' + iq C((q', q)) 



<5(k' + q + lq' C(q, q'))<J(p - q' - iq C((q', q)) 



(6.11) 



where the function iV(q', q) contains the product of two Jacobians and reads 
#(q,q') 



v/P(q) + P(q + q') ^(q') + P(q + q' 



4P(q + q') 



^(P(q)+P( q + q')) 2 -q' 2 ^/(P(q') + P(q + q')) 2 - q 2 



and the function C(q, q') is calculated as 

C(q',q) = l 



^/E(q)E(q>) 
E(q') - E(q' + q) 



E(q>) + E(q> + q) + VC^q 7 ) + P(q' + q)) 2 - q 2 



(6.12) 



(6.13) 



These permutation operators, which change the order of coupling, are essentially Racah coefficients for the Poincare 
group. In the nonrelativistic case the functions JV(q',q) and C(q',q) both reduce to the constant 1 and have 
the relatively compact form of the matrix elements of P given in e.g. [HI, [II]. Since both functions depend on 
magnitudes as well as angles between the momentum vectors, the 3D formulation is very appropriate for our relativistic 
calculations. In order to illustrate the momentum and angle dependence we display in Fig. Q] the function C(q',q) 
for a given value of |q'| = 0.65 GeV as function of |q| and several values of the angle y — q' • q. In general the values 
of C drop below 1 as q increases. The angle dependence is strongest for small q, where for q' ■ q — —1 the function is 
larger than 1. For the same fixed value of |q'| we display the function N(q' , q) in Fig. O Here we see a slowly varying 
dependence on the momentum |q| and a strong angle dependence. For small angles (y = 1) the function iV(q',q) is 
larger than 1, whereas for large angles (y = — 1) it is reduced from 1 by as much as 20%. 
In matrix form the Faddeev equation, Eq. (|5.6j) . has the form 

(6.14) 



(k,q||T||^,q ) = (k,q||T 1 P||^,q ) + (k,q|jT 1 P(z-Mo)- 1 T||^,qo). 



Since at this stage we are only carrying out first order calculations, we only need to consider the first term. Explicitly 
this is given as 

(k,q||T(W)||^,q ) = (k,q||T 1 (q,e)P|| V(I ,qo> 

= J dk'dq'dk"d q "(k, q ||T 1 ( q ,e)||k',q')(k',q'!|P||k",q")(k", q "|^, qo ) 

= J dk'dq'dk"dq" Ti(k, k', q; e)*(q' - q)<^(k")<5(q" - qo) 

xN(q', q")[5(k' - k(q", -q' - q"))<5(k" + k(q', -q' - q")) 
+S(k> + k(q", -q' - q"))*(p" - p(q', -q' - q"))] 

(6.15) 



N(q,qo) T s p, qo + -q C(q , q), q; e ) <p d [q+ -q C(q,q ) 



Here the invariant parametric energy e which enters the two-body t-matrix is given by 
consider bosons, the symmetrized two-body transition matrix T s is given by 

r s (k,k', q ;e) = T 1 (k,k', q ; £ )+T 1 (-k,k', q;£ ) 
= T 1 (k,k', q ; £ )+T 1 (k,-k', q;£ ). 



W — y m 2 + q 2 ■ Since we 



(6.16) 



This expression is the starting point for our numerical calculations of the transition amplitude in first order. The first 
step for an explicit calculation is the selection of independent variables. Since we ignore spin and iso-spin dependencies, 
the matrix clement (k, q||T||i^d; qo) is a scalar function of the variables k and q for a given projectile momentum qo. 
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Thus one needs 5 variables to uniquely specify the geometry of the three vectors k, q and qo. We follow here the 
procedure from Ref. [l§] and choose as variables 

k = l k |, 9 = |q|> x k = k • qo, x q = q ■ q , xf q = (q x q) ■ (q x k). (6.17) 

The last variable, , is the angle between the two normal vectors of the k-q -plane and the q-qo-plane, which are 
explicitly given by 

/ — - n x qo x k 

(q x k) = 



1 - (qo • k) 2 

(qo^q) = (6-18) 

V 1 - (qo ■ q) 2 

With these definitions of variables the expression for the transition amplitude as function of the 5 variables from 
Eq. (|6.17|) has the same form as its nonrelativistic counterpart, and we can apply the algorithm developed in Ref. [l2j 
for solving it without partial wave decomposition. The only additional effort is the evaluation of the functions 
C(go: Q, x q ) and N(qo, q, x q ). However, this is conceptually and numerically quite straightforward, since both functions 
depend only one angle, x q . 



VII. TWO-BODY TRANSITION OPERATOR AND RELATIVISTIC DYNAMICS 



The kinematic effects related to the use of the relativistic Racah coefficients have been described in the previous 
Section, ft is left now, to obtain the transition amplitude of the 2N subsystem embedded in the three-particle 
Hilbert space , T s (k, k', q; e), entering Eq. (|6.15[) . This is a fully off-shell amplitude depending in addition on the 
Poincare-Jacobi momentum q of the pair. The embedded 2N transition amplitude satisfies the Lippmann-Schwinger 
equation 



T 1 (k,k / ;q)=V 1 (k,k';q) 



dk" 



y 1 (k,k";q)T 1 (k",k';q) 



,2 _|_ p /2]2 



(7.1) 



t 2 + p" 2 ] 2 



where the interaction (|3.24[) can be expressed in the relativistic Jacobi momenta as Ref. [32 



V(q) = \ 2vm 2 + k 2 +v 



2VTO 2 



k 2 



q 2 . 



(7.2) 



For q = this expression reduces to the interaction v(k, k')<5(q— q'), which is the interaction in the two-nucleon mass 
operator. In the same limit, Eq. (|7.1[) reduces to the familiar Lippmann-Schwinger equation with relativistic kinetic 
energies. 

This matrix clement is constructed using the methods outlined in equations (|5.12[) and (|5. 14|) . First the matrix 
element of the right half-shell embedded t-operator is evaluated using the two-body half shell transition amplitude 
where the convention of [35| is employed 



(k|r 1 (q;z / )|k / ) = (k|nq)|k'(->) 

2(E k , +E k ) 



-Mk,k';2E k ,), 



y/4E%, + q 2 + v/4£ 2 + q- 

where the 2N transition amplitude t(k, k'; 2Ek>) is the solution of the half-shell Lippmann-Schwinger equation 

v(k,k")t(k",k';2E k ,) 



t(k,k';2E k ,)=v(k,k')+ / dk 



E k , - 2Vr 

This solution is used as input to equation (|5. 14)) which has the form 
(k|T 1 (q J z)|k / ) = (k|r 1 (q;z , |k / ) - 
dk" (klT^q^JIk" 



fc" 2 



(7.3) 



(7.4) 



v /4(?n 2 + k" 2 ) + q 2 



v /4(m 2 + k" 2 ) + q 2 



^"1^(^^)1^,(7.5) 
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Note that in this equation the unknown 



where Ti(z') is taken to be right half-shell with z' = y/ 4(m 2 + k' 2 ) + q 2 - 
matrix element is to the left of the kernel. 

We refer to Ti(k, k';q) := (k|Ti(q; z)|k') as the embedded 2N t-matrix and to V(k, k';q) as the embedded inter- 
action. Matrix elements of T\{z) can be alternatively calculated by inserting a complete set of eigenstates of the 2N 
mass operator myi = 2\/k 2 + m 2 + v, as has been carried out in Ref. [3^ | for a relativistic calculation of the triton 
binding energy using two-body s-waves. Similarly, this method of spectral decomposition can be used to directly 
calculate the matrix elements of the embedded two-body t-matrix, as has been done in another relativistic calculation 
of the triton binding energy (4p| . The general difficulty with this method of spectral decomposition is that it requires 
an integration over the cm. half-shell matrix elements t(k,k'; Ek>) in k', requiring the knowledge of those matrix 
elements for large values of k' , which can pose a challenge with respect to numerical accuracy. To our knowledge, this 
method has not yet led to a successful relativistic calculation of scattering observables. 

We use Eq. (|7.5p to explicitly construct the elements of the fully off-shell t-matrix, which enters the calculation of the 
three-body transition amplitude given in Eq. (|6 . 1 5[) . For every off-shell momentum k' the integral equation, Eq. (|7.5D , 
must be solved for each z. It is worthwhile to note that the k" integration in Eq. (|7.5[) only involves momenta of the 
half-shell t-matrices, but no energies. The momenta k and k' are fixed by requirements of the three-body calculation, 
and typically are not higher than 7 GeV. We tested that for converged results the k" integration has to go up to 
about 12 GeV. The singularities in the two denominators of Eq. (|7.5|) do not pose any problems and are handled with 
standard subtraction techniques. 

In order to obtain insight into the impact of the embedding for different values of q, we introduce approximations 
to the embedded interaction. First, we completely neglect q in the embedded interaction, which leads to 



V(k,k';q)^V (k,k';q)=v(k,k'). 
Furthermore, we want to test the leading order terms in a q/m and v/m expansion as suggested in Ref. [R 

V(k, k'; q) -> V 1 (k, k'; q) = w(k, k') (\ - -^j 

and 

y(k,k';q)^^ 2 (k,k';q) = ? ;(k,k) 
and explore their validity as function of projectile energy. 



(7.6) 



(7.7) 



8\/m 2 



k 2 y/m 2 



(7. 



VIII. CROSS SECTIONS FOR ELASTIC SCATTERING 

The calculation of the cross section for elastic scattering, Eq. (|6.5|) . requires the knowledge of the matrix element 
(ipd, q, qo\U\(pd, qo). Using the definition of the operator U, Eq. (|5.5|) . inserting a complete set of states and using for 
the matrix elements of the permutation operator the expression from Eq. (|6.11j) , we obtain 



(<Pd,q\\U\\<p d ,qD) = (<Pd,q\\P(z-Mo)+PT\\<p d ,qo) 

= J dk'd q 'dk"dq"(cp d , q|k', q')(k', q'||P||k", q")(k", q"|| (z - M) + T\\cp d , q ) 



= [W 



q 2 - W4(m 2 + (q + -q C(q, q )) 2 ) + qg 



x2V(q,q )<^d 
+2 J d 3 q" V(q,q> d 



Yd 



x (q+ -q" C(q,q"),q"||T||^,q ). 



q+ ^qo c(q, q ; 



q" + iqC(q",q) 



(8.1) 



In first order the transition amplitudes reads T = tP, thus the final expression for the transition amplitude for elastic 
scattering becomes 



(^,q||[/||^,q ) = ( W - x /m 2 + q 2 - \/4[m 2 + (q + -q C(q, q )) 2 ] + q 2 



1 
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x2N(q,q )ip d 
+2 J dq" N(q,q")cp d 



qo + ^<lC(qo,q) 



Yd 



q+ 2^0 C(q, q ; 



1 



q" + .qC(q",q; 



Yd 



q"+ 2qC(q",qo; 



(q + ^q" C(q, q")), (qo + C(q , q)); e 



(8.2) 



where e = W — ^/ m 2 — q" 2 . 

In the following we want to compare a non-relativistic first order calculation to a corresponding relativistic one. 
What is common in both calculations is the input two-body interaction. In the relativistic case it is transformed to 
be two-body scattering equivalent to the non-relativistic two-body calculation. Though we consider only spin-isospin 
independent interactions, we nevertheless can provide qualitative insights for various cross sections in three-body 
scattering in the intermediate energy regime, which we define as ranging from 200 McV to 1 GeV projectile energy. 
The focus of our investigation will be how kinematic and dynamic relativistic effects manifest themselves at different 
energies and for different scattering observablcs. 

As model interaction we choose a superposition of two Yukawa interactions of Malfhet-Tjon type [4l[ with pa- 
rameters chosen such that the potential supports a bound-state, the 'deuteron', at -2.23 MeV. The parameters are 
given in Ref. [l2j ■ With this interaction we solve the non-relativistic Faddeev equation in first order as a basis for all 
comparisons. Then we need to construct a phase equivalent relativistic two-body interaction. We use the procedure 
suggested by Kamada-Glockle [lH, and obtain a two-body interaction v(k, k') as Born term of a relativistic two-body 
Lippmann-Schwingcr equation. This two-body t-matrix, t(k, k';e) is the starting point for all calculations which will 
be presented in the following. In principle there are other methods to obtain a phase-shift equivalent relativistic poten- 
tial [20| . however in this work we want to focus on the relativistic effects visible in three-body scattering observablcs, 
and thus use only one fixed scheme. 

Following Ref. [l3l ]. as first assessment of the quality of different approximations for the embedded interaction we 
solve the relativistic 2N Schrodinger equation for the deuteron as function of the momentum |q|, which takes the form 



$«t(k) = — 1 = fdk' V(k,k';q)$ d (k'). (8.3) 

^W-^ 2 m +q 2 J 

Here nid is the rest mass of the deuteron. In Fig. [3] we show the deuteron binding energy e d calculated using the 
approximations of the embedded interaction given in Eqs. (|7.6[) . (|7.7p , and (|7.8|) . A correctly embedded interaction 
should of course not change Ed at all. We see that Ed based on the calculation using V starts to deviate already 
for very small q. The approximation V\ gives reasonable results up to q « 0.3 GeV, whereas V<z is good to about 
0.6 GeV. In the following we will see how far these simple estimates arc reflected in the calculation of various scattering 
observables. 

As a first observable we consider the total cross section for clastic scattering, er e ; , which is given in Table U for 
projectile kinetic energies from 10 MeV up to 1 GeV. Starting from the non-relativistic cross section, we successively 
implement different relativistic features to study them in detail. First we only change the phase space factor in the 
calculation (psf ) together with the relativistic transformation from laboratory to cm. frame, and only then implement 
the relativistic kinematics due to the Poincare-Jacobi coordinates (R-kin). The relativistic phase space factor alone 
has a large effect on the size of the total cross section, as was already observed in [13] ■ The kinematic effects of the 
Poincare-Jacobi coordinates have the opposite effect and lower the cross section. However, all kinematic effects taken 
together increase the total cross section by about 6% at 0.2 GeV and about 40% at 1 GeV. Introducing relativistic 
dynamic effects into the calculation changes this considerably . The full relativistic calculation (R) lowers the total 
cross section by about 2% at 0.2 GeV and about 6% at 1 GeV, so that in total the relativistic cross section is smaller 
than the nonrelativistic one. The approximation V2 of Eq. (|7.8p is very good in the energy regime considered, even 
at 1 GeV its result only deviates by about 2% from the full one. As suggested by the calculations of the deuteron 
binding energy, the approximation V\ of Eq. (|7.7p is still reasonable at 0.2 GeV, but after that starts to become worse. 

Next we consider the differential cross section for elastic scattering. In Fig. [4] we show the calculation for 0.2 GeV 
projectile kinetic energy. Since differences between the calculations disappear on a logarithmic scale, we also show 
the quantity 



\dflJ R 



NR 



> dQ ) NR 



(8.4) 
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expressed in percentage for the different approximations in the lower panel of Fig. [4j For the backward angles, 
9 > 135°, which correspond to higher momentum transfer, all relativistic effects increase the cross section. Here it 
can be clearly seen that indeed Vq is a bad approximation, whereas V\ and V 2 are of about the same quality. We 
also see that there is a small difference between the calculation based on V 2 and the full result. Similar findings, 
however without the full calculation, were presented in Ref. When going to higher projectile kinetic energies, 

we expect that the effects increase. This is indeed so, as shown in Fig. [5] for the differential cross section at 0.5 GeV 
projectile kinetic energy Here the second minimum in the cross section around 90° shows a shift towards larger 
angles once relativistic dynamics is included. This phenomenon has been seen and studied in some electron-deuteron 
scattering [43[ calculations. To study this shift in more detail we show in the lower panel of Fig. [5] a restricted angular 
range. Here we can see that the relativistic kinematics produces a shift of the minimum by a few degrees. The 
magnification shows that the approximations of the embedded interaction oscillate by a few degrees around the full 
solution. At the extreme backward angles, the relativistic cross section is larger than the non-relativistic one, as was 
the case at 0.2 GeV. In order to illuminate the details of the two minima of the cross section at 0.5 GeV even further, 
we show in Fig. [6] the two terms contributing to the operator U for elastic scattering separately. The curves labeled 
' 1st- ■IT correspond to the first term in Eq. (|8.2p or the operator P(z — Mo) in Eq. (|5.5j) . which contributes to the 
structure of the cross section at backward angles and depends only on the product of two deuteron wave functions 
evaluated at shifted momenta. Here the relativistic calculation is pushed slightly towards smaller angles indicating the 
effect of the functions C(q, qo). The second term in Eq. (|8.2[) . represented by the curves labeled l int-lT, contains in 
integral over the two-body t-matrix and a product of deuteron wave functions, and basically determines the structure 
of the cross section for angles up to about 100°. Here we see the shift of the minimum towards higher angles for the 
relativistic calculation. The interference of both terms in the calculation gives the final pattern as seen in Fig. [5] 

IX. CROSS SECTIONS FOR BREAKUP PROCESSES 

The calculation of the breakup cross section, Eq. (|6.4p . requires the knowledge of the matrix element 
(k, q||{/o||0d,qo)- Energy conservation requires that in Eq. (|6.4jl W/ = W, = W = y/ 4(m 2 + k 2 ) + q 2 + \J m 2 + q 2 . 
This gives a relation between the momenta k and q. In fact, for each given q the magnitude of k is fixed as 



k a = \k a \ = 5\/W 2 - 3m 2 - 2W v / m 2 + q 2 . This leads to 

CT ™ = ^q^ J dn p dn q dqq 2 ^^(m 2 + k 2 a )+q 2 | <p, q|| U \\<p d , q ) ? (9.1) 

We will consider here the cross sections for two different breakup processes, the inclusive breakup, where only one 
of the outgoing particles is detected, and the full or exclusive breakup. In order to obtain the differential cross section 
for inclusive breakup, one still needs to integrate over the solid angle of the undetected particle to arrive at invariant 
cross section: 

g^g _ (2^) 4 E(q )E d (q )E(q)q 



dQqdEq 3 ?oW j dQk^V^TkJj+q 2 |(k,q||C/ ||^,qo)r. (9.2) 

Here we changed from the variable dq to the more utilized dE q . The five- fold differential cross section for exclusive 
scattering, where both particles are detected, is given by 

d 5 <rF™ (2^) 4 E(q )E d (q )E(q) q p a 



dn k dfl q dE q 3 q \N 4 



^/4(m 2 + k 2 a )+q 2 |(k,q||[/ ||^,qo)r. (9.3) 



Next we need to explicitly evaluate the matrix element for breakup scattering, (k, q| | ?7o 1 1 7 Qo) ? with Uq given in 
Eq. (jST5j) 

(k,q||[/o||<Aj,q ) = (k ) q||T|| ¥ > d) qo) + (k,q||P 12 P 2 3T||^ ) qo) + (k ) q||Pi 3 P23T||^,qo). (9.4) 

The two terms containing the permutations can be calculated analytically, as we show explicitly for the second term 
using the expressions of Eqs. (|6.8p and (|6.10[) for the Poincare-Jacobi coordinates 

(k,q||P 12 P 23 T||^,q ) 
= J dqidq 2 dq 3 J dk'dq'(k, q|qi, q 2 , q 3 ) (qi, q 2 , q3||pL2P23||k', q') (k', q'||T|| tp d , q ) 
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= J dq 1 dq 2 dq 3 J dk'dq'<5(qi + q 2 + q 3 - P) n(q\; q 2 , q 3 ) n(q 2 ; q 3 , qi) 

xS(q - qi ) <5(k - k 23 ) 5(q' - q 2 ) 6(k' - k 31 )(k', q'\\T\\cp d , q ). (9.5) 
Taking particle 1 as spectator we can evaluate the momenta q^ explicitly as 

qi = q 

q 2 = k 



2E(k) I ^2EW) 



2 I „2 




q 3 = k ttt , +S(fc) . (9.6) 

2S(fc) ^(2^(fc)) 2 +9 2 + 2i?(fc) V V V j 

From this q' and k' can be obtained as q' = q 3 and k' = k 3 i by inserting the expressions of Eq. (|9.6j) into Eq. (|6.8p 
leading to 

(k,q||_Pi 2 P 23 T||<^,q ) = ra(qi;q 2 q 3 ) h(q 2 ; q 3 qi)(k 3 i, q 2 ||T||^ d , q ). (9.7) 
In first order we have T = tP, and an explicit evaluation leads to 

(k,q|| P 12 P 23 T y d ,qo) = 

-( \«t \ ati \rr (i , C '(q0'q2) \ ( , C(q 2 ,q ) \ 

?i(qi;q 2 q 3 ) n(q 2 ;q 3 qi) A^(q 2 ,q ) T s k 3 i,q H q 2 ;e <p d q 2 H q .(9.8) 



The functions A(q 2 ,qo) and C(qo,q 2 ) are defined in Eqs. (|6.12[) and (|6.13[) . The last term in Eq. 
(k, q||Pi 3 P 23 T||(^ ( i, qo), is calculated analogously. Having calculated the matrix element of Uq, Eq. (|9.4[) . we can 
obtain the differential cross section for inclusive as well as exclusive breakup scattering. The expressions for the 
invariant cross sections in the laboratory variables are derived in Appendix [S] 

First, we consider inclusive breakup scattering and compare the cross sections for a non-relativistic first-order 
calculation in the two-body t-opcrator with the corresponding relativistic one. One can expect that the evaluation of 
the delta function in the cross section, Eq. (|4.20| will have a substantial effect on breakup cross sections, since it fixes 
the relation between the magnitudes of the vectors k and q. This in turn determines the maximum allowed kinetic 
energy the ejected particle is allowed to have as function of the emission angle. To get a global impression of those 
differences Fig. [7] shows a contour plot of the differential cross section for inclusive breakup scattering as function 
of the kinetic energy and the emission angle of the ejected particle for the non-relativistic and the fully relativistic 
calculation. The figure shows that for each angle the maximum allowed kinetic energy of the ejectile is shifted in 
the relativistic calculation towards smaller values compared to the non-relativistic calculation. Specifically, one can 
expect a shift of the quasi-free scattering (QFS) peak usually studied in inclusive breakup scattering experiments. In 
Figs. [8] and [9] we present specific cuts at different constant angles to study details of the calculation. The upper panel 
of Fig.[5]shows the entire energy range of the ejectile at emission angle 6\ = 24° in a logarithmic scale, while the lower 
two panels give a close-up of both peaks on a linear scale. The QFS peak at the large ejectile energy exhibits clearly 
a shift towards a slightly lower energy compared with the peak position calculated non-relativistically. At this angle, 
relativistic kinematics given by the phase-space factor and the Poincare-Jacobi coordinates and indicated by the line 
labeled 'R-kin' results in a peak height, which is almost double that of the full relativistic calculation shown as solid 
line labeled 'R'. For breakup scattering we also study the different approximations to the full calculation as introduced 
in Section IV. In the QFS peak, which is defined by the condition that one of the particles is at rest, even the crudest 
approximation Vq, Eq. (|7.6[) . is very close to the full calculation, the approximations V\ and V 2 are indistinguishable. 
This is not surprising, since having one particle at rest means that the remaining two are almost in their own cm. 
frame, thus 'boost' effects should be extremely small. Note that we work here in the total cm. frame. It is quite 
illuminating to consider the second peak at very small ejectile kinetic energies. Since the energy of the particle is 
very small, it should become essentially non-relativistic. This is indeed the case, and the full relativistic calculations 
almost coincides with the non-relativistic one. It is worthwhile to note that neither relativistic kinematics alone nor 
the approximation Vq, which neglects the dependence of the embedded interaction on the pair- momentum is close to 
the non-relativistic and full relativistic calculations. However, an approximate consideration of this dependence as 
given by V\ or V 2 seems to suffice. We found that this behavior is similar for low energy ejcctilcs, independent if the 
energy of the projectile is 0.5 GeV or 1 GeV. In Fig. [5] we show the QFS peak for two different angles in order to 
convey that the increase or decrease of the height of the peak depends on the emission angle under consideration. In 
Fig.[TO]we show the QFS peak calculated for a projectile energy of 0.495 GeV and emission angles of 18° and 24°, since 
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there is experimental information available for one of the angles. Here we see that the relativistic calculation puts the 
peak at a position consistent with the data. Since we only carry out a first order calculation with a model potential, 
we are not surprised that the height of the peak is not described. A similar observation concerning the peak position 
was already made in Ref. [44|, where a first order calculation with two realistic NN interactions was carried out. To 
give an indication how the position of the QFS peak shifts with increasing projectile energy we show in Fig. [IT] the 
inclusive cross section for an emission angle of 24° for projectile energies 0.8 GeV and 1.0 GeV. We see again that in 
the QFS regime the approximations V% and Vi are essentially indistinguishable from the full calculation. Considering 
only effects of relativistic kinematics results in a peak height double as large as the full calculation, indicating that 
dynamic effects are very important at those high projectile energies. 

When considering exclusive breakup one faces many possible configurations that could be considered. Since we are 
carrying out a model study, we only want to show three specific configurations at two selected energies, 0.5 GeV and 
1 GeV in Figs. H"2l and [T51 In all three configurations the angle cfikq between the projections of the vectors q and k 
into the plane perpendicular to the beam direction qo is zero. Naive expectation is that when considering scattering 
in first order in the i-operator the probability that one of the particles is scattered along the beam is large. That 
is indeed the case as shown in the upper panels of Figs. [T2l and IT3| depicting a so called collinear configuration in 
which the angle between the vector q and the beam direction qo is zero (q ■ qo = x q = 1) and the angle between 
the vector p and the beam direction is 90° (k ■ qo = Xk — 0). Once the collinear condition is no longer fulfilled, 
the cross section becomes considerably smaller, as can be seen in the middle and lower panels of Figs. [T2] and [T3"l 
For the middle panels the angles are given by x q = V3/2 and x p = —0.5, for the lower panels they are x q = —0.25 
and Xk = —0.9. All configurations in Figs. [T2l and [TBI show considerable difference between the nonrelativistic and 
relativistic calculations. At 1 GeV we were specifically looking for possible configurations where the approximations 
V\ and V2 of Eqs. (|7.7p and (|7.8[) are not close to the full relativistic calculation any longer. Considering the two-body 
binding energy displayed in Fig. [31 on should expect that the approximation V2 can exhibit deviations from the full 
result at 1 GeV. One such configuration is shown in the lower panel of Fig. [T31 where there is a big difference between 
the calculation with V\ and the full result and a discernible difference between the calculation based on V2 and the 
full result. However, we note, that in the majority of configurations where the cross section is still relatively large, 
V2 is still a very good approximation at 1 GeV. Looking at the same configuration at 0.5 GeV, V2, and even V\ are 
extremely good approximations to the full result. 



We investigated relativistic three- nucleon scattering with spinless interactions in the framework of Poincare invariant 
quantum mechanics. Since that framework is not widely used in the nuclear physics few-body community we thought 
it adequate to discuss the formulation of that scheme in some detail, as well as the formulation of scattering theory 
in this framework. The main points are the construction of unitary irreducible representations of the Poincare group, 
both for noninteracting particles as well as for interacting ones. The Poincare interacting dynamics is constructed 
by adding an interaction to the noninteracting mass operator which commutes with and is independent of the total 
linear momentum and the z-component of the total spin. 

In this work we do not use partial waves but rather internal vector variables. This leads to what we call Poincare- 
Jacobi momenta. They Wigner rotate under kinematic Lorentz transformations of the underlying single particle 
momenta. In the interacting three-body mass operator the two-body interactions are embedded in the three-particle 
Hilbert space and are given as [32j 



This expression shows explicitly the dependence on the total momentum of the two-body system. For the sake of 
completeness we also discuss the multichannel scattering theory in that relativistic framework and established the 
manifestly invariant form of the differential cross section. 

The application to three-body scattering is based on the Faddeev scheme, which is reformulated relativistically 
working with various types of mass operators. The usage of the Poincare-Jacobi momenta leads to algebraic modifica- 
tions of corresponding standard nonrelativistic expressions, like e.g. the momentum representation of the permutation 
operator, Jacobians for the transitions between individual Jacobi momenta. Due to the dependence on the total mo- 
mentum of the two-body interaction embedded in the three-body system, the two-body off-shell t-opcrator entering 
the Faddeev equation acquires additional momentum dependence beyond the usual energy shift which is characteristic 
in nonrelativistic calculations. This two-body t-operator is then evaluated by expressing it exactly in terms of the 
solution of half-shell Lippmann-Schwingcr equations for a given two-body force in its cm. frame. We also solve 



X. SUMMARY AND CONCLUSIONS 




(10.1) 
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the relativistic Lippmann-Schwingcr equation for three different momentum dependent two-body forces, which are 
approximations to the relativistic embedded interaction. 

In order to compare a nonrclativistic calculation to a relativistic one, scattering equivalent two-body forces in the 
relativistic and nonrelativistic formulations have to be used. In this work we follow the KG prescription to arrive 
at scattering equivalent two-body forces. There are different schemes, and a detailed study on differences between 
those schemes will be subject of a forthcoming work. We also restrict ourselves to a first order calculation in the 
two-body t-operator, which however, already exhibits most of the new relativistic ingredients , both kincmatically 
and dynamically. The two-body force was chosen as a superposition of two Yukawas of Malflict-Tjon type supporting 
a bound state, the deuteron, at -2.23 MeV. We calculate three-body scattering observablcs in the intermediate energy 
regime, which we take to range from 0.2 GeV to 1 GeV. Those observables are cross sections for elastic as well as 
breakup scattering, namely inclusive and exclusive scattering. Not surprisingly we find that the difference between 
nonrelativistic and relativistic calculations increase with increasing energy. This is specifically apparent when looking 
at the positions of minima in the differential cross section as well positions of QFS peaks in inclusive scattering. When 
studying various approximations to the relativistic embedded interaction, we find that if the approximations contain 
terms proportional to the first order in a p/m and q/m expansion, the approximation captures the features of the exact 
relativistic calculation very well. Only at 1 GeV we start to find discernible deviations in selected configurations for 
exclusive scattering. Our results clearly indicate an interesting interplay of kincmatically and dynamically relativistic 
effects, which as expected, increase with energy. For example, the total cross section is increased by kinematical 
effects, whereas the dynamical effects resulting from the q-dependence of the embedded two-body force decrease it. 
This tendency is also seen in the height of the QFS peak where in most cases the full calculation is lower than the one 
allowing only for relativistic kinematic effects. Thus, considering only relativistic kinematic effects leads in general to 
an over prediction of cross sections, which becomes more dramatic the higher the projectile energy is. 

The Poincare invariant relativistic framework is formally close to the nonrelativistic one and therefore standard 
nonrclativistic formulations, in our case the Faddeev scheme, can be used with proper modifications. The present re- 
striction to a first order calculation in the t-operator will soon be replaced by a complete solution of the corresponding 
relativistic Faddeev equation following [l2j|, where the nonrelativistic Faddeev equation has been successfully solved 
employing vector variables. The application to the realistic world of pd scattering at the energies up to 1 GeV consid- 
ered in this study requires of course two-and three-nucleon forces high above the pion production threshold. Though 
they are not yet available, they can also be included in the framework of Poincare invariant quantum mechanics. 
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APPENDIX A: INVARIANCE OF S AND RELATION TO T 



The expression (|4.20[) and (|4.2ip for the differential cross section can be rewritten in a manifestly invariant form. 
We write them as a product of an invariant phase space factor, an invariant factor that includes the relative speed, 
and an invariant scattering amplitude. 

To identify and establish the invariance of the invariant scattering amplitude note that the scattering operator S is 
Poincare invariant: 

U {A,a)S = U {A,a)fll(H,H )n_(H,H ) = fll(H,H )U(A,a)Q-(H,H ) 

= ftl(H,H )Q.{H,H )U (A,a) = SU (A,a). (Al) 

The Poincare invariance of the S operator above is a consequence of the intertwining relations for the wave operators 

U(A,a)n±(H,H )=n ± (H,H )U (A,a) (A2) 

To show the intertwining property of the wave operators first note that the invariance principle gives the identity 

n±(H, H ) = fi±(M, M ). (A3) 

The mass operator intertwines by the standard intertwining properties of wave operators. For our choice of irreducible 
basis the intertwining of the full Poincare group follows because all of the generators can be expressed as functions 
of the mass operator and a common set of kinematic operators, {K, j z , j x , j 2 , — iV/c}, that commute with the wave 
operators. 

The covariancc of the S matrix elements follows from the Poincare invariance of the S operator if the matrix 
elements of S are computed in a basis with a covariant normalization. 

The S-matrix elements can be evaluated in the channel mass eigenstates. After some algebra one obtains: 
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(A4) 



where M a \a) = \N a \a) and Mp\fj) = W^|a) and W := ^(W Q + Mg) is the average invariant mass eigenvalue of the 
initial and final asymptotic states. In deriving (|A4[) the two strong limits in (|4.6p are replaced a single weak limit. 
Equation (|A4[) is interpreted as the kernel of an integral operator. S'-matrix elements are obtained by integrating the 
sharp eigenstates in Eq. (|A4|) over normalizablc functions of the energy and other continuous variables. To simply 
this expression define the residual interactions V a and V 13 by: 



V a := M - M a ; V = M - Mp, 

where 

V a \a) = (M - W„)|o); V \/3) = (M - W fl )|/3). 
The resolvent operators of the mass operator and the channel mass operator, 
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arc related by the second resolvent relations [37^ 

G(z) ~ G a (z) - G a (z)V a G(z) = G(z)V a G a (z). 



(A5) 
(A6) 

(A7) 

(A8) 



It is now possible to evaluate the limit as e — » 0. It is important to remember that this is the kernel of an integral 
operator. 
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The first term in square brackets is unity when the initial and final mass eigenvalues are identical, and zero otherwise; 
however, the limit in the bracket is a Kronecker delta and not a Dirac delta function. For a ^ (3, (/3(W')|a(W)) are 
Lebesgue measurable in W for fixed W, so there is no contribution from the first term in Eq. (|A4[) . For the case that 
W Q = \Np, we have (/3(W')|a(W)) oc S(W - W). The matrix element vanishes by orthogonality unless = W Q , but 
then the coefficient is unity. Thus, the first term in (|A4[) is (/3\a) if the initial and final channels are the same, but 
zero otherwise. The matrix elements also vanish by orthogonality for two different channels governed by the same 
asymptotic mass operator with the same invariant mass. The first term in (|A4[) therefore includes a channel delta 
function. 

For the second term, the quantity in square brackets becomes —2iri5(\Np — W Q ), which leads to the relation 

(p\S\a) = (a\(3) - 2iriS(\N fj - \N a )(f3\T^ a (\N a + iO + )\a), (A9) 

where 

T 0a (z) = V a + V p G{z)V> 3 . (A10) 

and (a is zero if the initial and final channels are different and is the overlap of the initial and final states if the 
initial and final channels are the same. Equation (|A9[) is exactly eq. (|4.13[) . 

With our choice of irreducible basis the residual interactions and the resolvent commute with the total linear 
momentum operator, and if the sharp channel states |a) and \(5) are simultaneous eigenstates of the appropriate 
partition mass operator and the linear momentum, then a three-momentum conserving delta function can be factored 
out of the T-matrix element: 

(/3|T^(W Q + iO+)\a) = S 3 (P - P a )(p\\T^(\N a + iO+)\\a). (All) 

When combined with the three-momentum conserving delta function the invariant mass delta function can be replaced 
an energy conserving delta function 
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The 5*-matrix elements can be expressed in terms of the reduced channel transition operators as follows: 

(j3\S\a) = (a\p)6 Pa - i(2n)6 4 (P - P a )^L(/3\\T^ a (\N a + iO+)\\a) (A13) 

In this expression the S operator is invariant while the single particle asymptotic states have a non-covariant normal- 
ization. 

To extract the standard expression for the invariant amplitude the single particle states are replaced by states with 
the covariant normalization used in the particle data book [36j |: 

|p, M ) — > \p,(,) cov = \p,^)^2EZ(^) 3/2 - (A14) 

The resulting expression 

- i(2ir)S 4 (P l3 - P a )^cov(0\\T^(\N a + iO + )\\a) cov (A15) 

is invariant (up to spin transformation properties). Since the four dimensional delta function is invariant, the factor 
multiplying the delta function is also invariant (up to spin transformation properties). This means that 

cov (a\\M al3 \\(3) cov := ^—^ cov ((3\\T^ a (\N a + iQ + )\\a) cov (A16) 

is a Lorentz covariant amplitude. The factor of l/(27r) 3 is chosen to agree with the normalization convention used in 
the particle data book [36| . 

The differential cross section becomes 

(27r) 4 2 
d<7 = — — r \ cov (p l7 - ■ ■ ,p n ,\\M al3 \\p b ,p t ) cov \ 

4:E mt (pt)E mb (p b )v bt 1 1 

■"*(?».-*-*) 5^?- (A17) 
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The identity 



Vbt = v^gZ^M (A18) 

can be used to get an invariant expression for the relative speed between the projectile and target and 



d$n(Pb + h;pi, ■■ - ,Pn) = S 4 ^2 



is the standard Lorentz invariant phase space factor. Inserting these covariant expressions in the definition of the 
differential cross section gives the standard formula for the invariant cross section 

(27r) 4 2 

do- = \cov(pi, ■ ■ ■ ,Pn, \\M a/3 \\p b ,pt)c v\ d$ n (p b +p 4 ;Pi,- • - ,Pn)- (A20) 



4\/(pt -Pb) 2 - 



Because of the unitarity of the Wigner rotations and the covariance of | ccw (pi, • • • ,Pm \\M a ^\\pb, ptjcov | this becomes 
an invariant if the initial spins are averaged and the final spins are summed. 

This manifestly invariant formula for the cross section is identical to (|4.20j) and (|4.21| : in this form the invariant 
cross section can be evaluated in any frame. The index t refers to the target, which is in our case the deuteron. 



APPENDIX B: BREAKUP CROSS SECTION IN THE LABORATORY FRAME VARIABLES 



The total breakup cross section is Lorentz invariant. The expression for the differential cross sections (|4.20|) is given 
in terms of single particle variables, while the solutions of the Faddeev equations give transition matrix elements as 
functions of the Poincare-Jacobi momenta defined in Section II. To compute the total cross section it is useful to 
work in a single representation. Since the single particle momenta are directly related to measured parameters of 
the differential cross section, we change variable in the transition amplitudes from Poincare Jacobi moment to single 
particle momenta. In this section the single particle variables are computed in the laboratory frame. 

The relation between the product of single particle basis states and states expressed in terms of the Poincare-Jacobi 
momenta are given by 



|Pn,Pd; 



Ipi,P2,p3; 



9(qo,P,) 



d{Pn,Pd) 

9(k, P2 + 



1/2 



P3) 



|qo,P* 



1/2 



<9(P2,P3) 



3(q,P, 



<9(pi,p 2 + P3) 



1/2 



P/,k,q) 



(Bl) 



The Jacobians in these transformations are 



<9(qo,P,: 



d(p n ,Pd) 
<9(k, p 2 + p 3 ) 



E(g )E d (g )E i 
E{ Pn )E d {p d )W 



<9(P2,P3) 

0(q,P/) 



Defining 



<9(pi,P2 + P3) 



Wk,q) 



+ E( P3 ))2 - (p 2 + p 3 )~ (E( P2 ) + E(p 3 )) 
AE{p 2 )E{p 3 ) 

y/(E(jp2) + E{p 3 )f - (p 2 + paF+q* E(g)E f 
Efa.) (E( P2 ) + E{p 3 )) W 



9(qo,P») 

E? 



d(p n ,Pd 

the total cross section for breakup scattering becomes 

lab 1 (2VT) 4 



<9(k, p 2 + p 3 ) 



9(P2,P3) 



%P/) 



9(Pl.P2 + P3) 



(B2) 



(B3) 



f d Pl dp 2 dp 3 6(E f - Ei) S 3 (P f - Pi)T(W,p,q) \(p,q\\Uo\\<Pd,<is>) 

5— / dpidp2<5(E / -E J )r(iy,p,q)|<p,q||[/ ||^,qo>| 2 , 
>5 Vbt J 



(B4) 
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where we used P = Pj = P/. The S function in the energy can be eliminated by a variable change. The total energy 
of the system is E = E, = E/ = \Jm 2 + p\ + ^m 2 + p\ + \/m 2 + (P - p x - p 2 ) 2 . With pi = |pi| and p 2 = P2/IP2I 
this gives 



<7E 



P2 - (P - Pi) • P2 P 2 (E - - B(pa)(P - pi) • p 2 



dp 2 B(pa) E - E(p±) — Eipz) 
Since dpi = E(j>\) jp\dE\, and v b t =p n /E(p n ), Eq. (|F34|) becomes 



E(p2)E(p 3 ) 



Jab _ (2tt) 4 £;(p„) 



'br 



dnxdn 2 dE 



PlP 2 2 E{ Pl )E(p 2 )E(p 3 ) 



1 p 2 (E~E( Pl ))-E(p 2 )(P-p 1 )-p 2 



3 Pn 

Inserting the explicit expression of Eqs. (|B2[) and (|B3[) we obtain 

(2tt) 4 E(q )E d (q ) 



r(W > k,q)|{k ) q||Ub||^,qo> 



(B5) 



(B6) 



Jab 
'br 



2p n m d 



J dfl 1 dfl 2 dE 1 E{q)\/m 2 + kV4(m 2 + k 2 ) + q 2 



PiP 2 



j>a(E-j;(pi))-£7(pa)(P-pi)-pa 



r(W,k,q)|(k,q||f/ ||^,q ) 



(B7) 



This gives the total invariant cross section as a five dimensional integral. We have expressed it as a function of the 
incident laboratory momenta. 

It follows that the five-fold differential cross section for exclusive breakup scattering 



d 5 rr lab 



dttidn 2 dEi 



(2tt) 4 E(q )E d (q ) )A /^^^ 4(m2+k 2 )+ 2 

3 2p n m d 



> 2 (E - E{ Pl )) - E(p 2 )(P - pi) • p 2 



r(W,k,q)|(k,q||C/ ||^,qo) 



(B8) 



In inclusive breakup scattering only one of the outgoing particles is detected. Thus the cross section still contains 
an integration over the coordinates of the undetected particle. In order to calculate this cross section, it is convenient 
to start again from Eq. (|B4[) . However, since we need to integrate over the coordinates of one of the particles, we pick 
without loss of generality particle 1 as spectator and use as coordinates 



lab 
'br 



where we define 



1 (2tt) 4 



3 I'M 
1(2704 

3 v bt 



dpidkdp 
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<9(p 2 ,P3) 



<9(k,p 2 + p 3 ) 



6(E f - EO^Py - P 4 )r(W,k,q) Kk.qUC/oll^qo) 



dp 1 dfa5(E / -E i )r / (W > k > q)|(k,q||^ ||v«i,qD> 



<9(p 2 ,P3) 



<9(k,p 2 + p 3 ) 



Wk,q). 



(B9) 



(BIO) 



Since <5(E/ — E^) = -^S(\N ' j — Wj) and rfW jdk = 4k / '4(m 2 + k 2 ) + q 2 , the integration over k is eliminated leading 



to 



„lab 
a br 



3 V bt 



dn 1 dE 1 fe 1 E(fci)r / (w,k,q) j dOfcKk.qllc/oll^qo) 



Insert Eq. (|B10|) gives the explicit expression for the inclusive breakup scattering cross section 



iQb _ (2tt) 4 E(qo)E d ( qo ) 

a br ~ n a 

3 4p n m d 
and the differential cross section 



dflidE 



Pl kE{q)(4(m 2 + k 2 ) 



(P-Pi) 2 



d 3 a l b f _ (2tt) 4 E(q )E d (q ) p lP E{q){A(m 2 + k 2 ) + q 2 ) 
dVLidEx 3 Ap n m d ^4(m 2 + k 2 ) + (P- Pl ) 2 



dn fc |<k, q||C7b||v»«l, qto) 



dfl k \(k,q\\Uo\\<p d ,q.o) 



(Bll) 



(B12) 



(B13) 
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E lab [GeV] 


NR [mb] 


psf [mb] 


R-kin [mb] 


V [mb] 


Vi [mb] 


V 2 [mb] 


R [mb] 


0.01 


100027.1 


100766.0 


100605.6 


100363.2 


99288.5 


99394.2 


99276.9 


0.1 


398.5 


445.0 


443.9 


418.9 


397.5 


400.1 


399.2 


0.2 


167.5 


185.7 


184.5 


173.4 


163.2 


164.5 


164.1 


0.5 


67.6 


83.2 


81.8 


73.5 


63.9 


65.4 


65.4 


0.8 


42.9 


58.6 


57.1 


48.3 


38.7 


40.4 


40.8 


1.0 


34.2 


49.7 


48.1 


39.6 


29.7 


31.5 


32.3 



TABLE I: The total cm. cross section a for elastic scattering calculated from a Malfiiet-Tjon type potential. The nonrelativistic 
total cross section is given in the 2nd column, labeled NR, the relativistic one is given in the last column, labeled R. The other 
columns give the total cross section for elastic scattering when different relativistic features are successively implemented: psf 
shows the effect of the relativistic phase space factor, R-kin adds the relativistic kinematic effects resulting from the permutation 
operator, for Vb the relativistic two-body LS equation is solved with a cm. interaction, and Vi (i=l,2) denote the approximations 
of the embedded interaction given in Eqs. (|7.7[1 and (|7.8[) . 
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FIG. 1: (Color online) The momentum and angle dependence of the function C(q',q) from Eq. (|6.13|) at fixed momentum 
q' = 0.65 GeV. 



q' = 0.65 GeV 



q 1 




FIG. 2: (Color online) The momentum and angle dependence of the function iV(q',q) from Eq. (|6,12[) at fixed momentum 
q' = 0.65 GeV. 
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FIG. 3: (Color online) The deuteron binding in energy calculated with the embedded interaction ^(k, p; q) as function of 
q. The solid line labeled 'R' represents the binding energy of -2.23 MeV which is independent of q, when the full embedded 
interaction is employed. The dotted line is obtained if q is set to zero in the embedded interaction. The dotted, dash-dotted 
and dashed lines show the approximations to the embedded interaction Vo, Vi, and V2 as given in Eqs. (|7.6|l . (|7.7|l . and (|7.8|l . 
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FIG. 4: (Color online) The differential cross section for elastic scattering at 0.2 GeV projectile kinetic energy as function of the 
laboratory scattering angle (upper panel). The solid line represents the fully relativistic calculation. The lower panel shows 
the relative deviation A with respect to the corresponding non-relativistic calculation. The long-dashed curve labeled 'R-kin' 
represents a calculation in which only relativistic kinematic effects are incorporated. The dotted, dash-dotted, and dashed 
curves show the approximations to the embedded interaction Vo, Vi, and V2 as given in Eqs. (|7.6|l . (|7.7[) , and (17.81) . 
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FIG. 5: (Color online) The differential cross section for elastic scattering at 0.5 GeV projectile kinetic energy as function of 
the laboratory scattering angle. The double-dotted curve labeled 'NR' represents the non-relativistic calculation, and the solid 
curve labeled 'R' the corresponding fully relativistic one. The long-dashed curve labeled 'R-kin' represents a calculation in 
which only relativistic kinematic effects are incorporated. The dotted, dash-dotted, and dashed curves show the approximations 
to the embedded interaction Vo, Vi, and Vi as given in Eqs. (|7.6|l . (|7.7|l . and (|7.8|l . 
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FIG. 6: (Color online) The differential cross section for elastic scattering at 0.5 GeV projectile kinetic energy as function of 
the laboratory scattering angle. The double-dotted curve labeled 'NR' represents the non-relativistic calculation, and the solid 
curve labeled 'R' the corresponding fully relativistic one. The two other sets of curves show the contributions from the two 
different terms contributing to the transition operator U in the non-relativistic (NR) and relativistic (R) calculation. See text 
for further discussion. 
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E lab = 500 MeV 
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FIG. 7: (Color online) The inclusive cross section at 0.5 GeV laboratory projectile kinetic energy as function of the energy 
E of the emitted particle and its emission angle 9. The angles above the thin dashed line for larger energies and angles 
are kinematically forbidden. The thin lines represent the contours of the non-relativistic calculation, whereas the thick lines 
represent the corresponding contours of the relativistic calculation. 
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FIG. 8: (Color online) The inclusive cross section at 0.5 GeV laboratory projectile kinetic energy as function of the energy E of 
the emitted particle and a 24° emission angle. The upper panel displays the entire energy range of the emitted particle, whereas 
the two lower panels show only the low and the high energies in a linear scale. The solid line represents the fully relativistic 
calculation, the double-dotted line the corresponding non-relativistic one. The long-dashed curve labeled 'R-kin' represents a 
calculation in which only relativistic kinematic effects are incorporated. The dotted, dash-dotted, and dashed curves show the 
approximations to the embedded interaction Vo, Vi, and V2 as given in Eqs. (|7.6[1 . (|7.7p . and (|7.8p . 




FIG. 9: (Color online) Same as Fig. [8] but for different fixed angles of the emitted particle. 
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E 1qK = 0.495 GeV 
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FIG. 10: (Color online) The inclusive cross section at 0.495 GeV laboratory projectile kinetic energy as function of the energy 
E of the emitted particle and fixed emission angles of 18° and 24° degrees. The solid line represents the fully relativistic 
calculation and the double-dotted line the non-relativistic one. The long-dashed curve labeled 'R-kin' represents a calculation 
in which only relativistic kinematic effects are incorporated. The data are from Ref. [45[. 
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FIG. 11: (Color online) The inclusive cross section at 0.8 GeV (upper panel) and 1 GeV (lower panel) laboratory projectile 
kinetic energy as function of the energy E of the emitted particle and the fixed emission angle of 24° degrees. The notation of 
the curves is the same as in Fig. [9] 
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E lab = 0.5 GeV 
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FIG. 12: (Color online) The exclusive breakup cross section at 0.5 GeV projectile kinetic energy as function of the ejected 
particle kinetic energy for three different configurations defined in the three-body cm. frame. For all configurations the angle 
4> pq is 0°. For the upper panel the cos of the angle between q and the beam qo is x q — 1, i.e. the scattering occurs along the 
beam line, in the middle panel x q — V3/2 and in the lower panel x q = —0.25. The cos of the angle between p and qo is in the 
upper panel x p = 0, in the middle panel x v = —0.5 and in the lower panel x v = —0.9. 
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FIG. 13: (Color online) Same as Fig ll2l but for projectile energy 1.0 GeV. 



